All-optical switching, bistability, and slow-light transmission 
in photonic crystal waveguide-resonator structures 
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We analyze the resonant linear and nonlinear transmission through a photonic crystal waveg- 
uide side-coupled to a Kerr-nonlinear photonic crystal resonator. Firstly, we extend the standard 
coupled-mode theory analysis to photonic crystal structures and obtain explicit analytical expres- 
sions for the bistability thresholds and transmission coefficients which provide the basis for a detailed 
understanding of the possibilities associated with these structures. Next, we discuss limitations of 
standard coupled-mode theory and present an alternative analytical approach based on the effective 
discrete equations derived using a Green's function method. We find that the discrete nature of the 
photonic crystal waveguides allows a novel, geometry-driven enhancement of nonlinear effects by 
shifting the resonator location relative to the waveguide, thus providing an additional control of res- 
onant waveguide transmission and Fano resonances. We further demonstrate that this enhancement 
may result in the lowering of the bistability threshold and switching power of nonlinear devices by 
several orders of magnitude. Finally, we show that employing such enhancements is of paramount 
importance for the design of all-optical devices based on slow-light photonic crystal waveguides. 

PACS numbers: 42.65.Pc; 42.70.Qs; 42.65. Hw; 42.79.Ta 



I. INTRODUCTION 

It is believed that future integrated photonic circuits 
for ultrafast all-optical signal processing require different 
types of nonlinear functional elements such as switches, 
memory and logic devices. Therefore, both novel physics 
and novel designs of such all-optical devices have at- 
tracted significant research efforts during the last two 
decades, and most of these studies utilize the concepts of 
optical switching and bistability [1]. 

One of the simplest bistable optical devices which can 
find applications in photonic integrated circuits is a two- 
port device which is connected to other parts of a circuit 
by one input and one output waveguide. Its transmis- 
sion properties depend on the intensity of light sent to 
the input waveguide. Two basic realizations of such a 
device can be provided by either direct or side-coupling 
between the input and output waveguides to an optical 
resonator. In the first case, we obtain a system with res- 
onant transmission in a narrow frequency range, while in 
the second case, we obtain a system with resonant reflec- 
tion. Both systems may exhibit optical bistability when 
the resonator is made of a Kerr nonlinear material. The 
resonant two-port systems of the first type, with direct- 
coupled resonator, can be realized in one-dimensional 
systems, and they have been studied in great details in 
the context of different applications. In contrast, the res- 
onant systems of the second type, with side-coupled res- 
onators, can only be realized in higher-dimensional struc- 
tures, and their functionalities are not yet completely un- 
derstood. 

Our goal in this paper is to study in detail the sec- 



ond class of resonant systems based on straight optical 
waveguides side-coupled to resonators as shown in Fig. 1. 
Moreover, we assume that the waveguide and resonator 
are created in two- or three-dimensional photonic crystal 
(PhC) [2], Due to a periodic modulation of the refractive 
index of PhCs, such structures may possess complete pho- 
tonic band gaps, i.e. regions of optical frequencies where 
PhCs act as ideal optical insulators. Embedding carefully 
designed cavities into PhCs, one can create ultra-compact 
photonic crystal devices which are very promising for ap- 
plications in photonic integrated circuits. As an illustra- 
tion, side-coupled waveguide-resonator systems created 
in PhCs through arrays of cavities are schematically de- 
picted in Fig. 1(b) and Fig. 1(c). 

Practical applications of such PhC devices are becom- 
ing a reality due to the recent experimental success in 
realizing both linear and nonlinear light transmission in 
two-dimensional PhC slab structures where a lattice of 
cylindrical pores is etched into a planar waveguide. In 
particular, Noda's group have realized coupling of a PhC 
waveguide to a leaky resonator mode consisting of a 
defect pore of slightly increased radius [3-6]; Smith et 
al. demonstrated coupling of a three-line PhC waveg- 
uide with a large-area hexagonal resonator [7]; Seassal 
et al. have investigated the mutual coupling of a PhC 
waveguide with a rectangular microresonator [8]; No- 
tomi et al. [9] and Barclay et al. [10] have observed 
all-optical bistability in direct-coupled PhC waveguide- 
resonator systems. 

Photonic-crystal based devices offer two major advan- 
tages over corresponding ridge-waveguide systems: (i) 
the PhC waveguides may have very low group veloci- 
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FIG. 1: (Color online) Three types of the geometries of a 
straight photonic-crystal waveguide side-coupled to a non- 
linear optical resonator, A a . Standard coupled- mode the- 
ory is based on the geometry (a) which does not account 
for discreteness-induced effects in the photonic-crystal waveg- 
uides. For instance, light transmission and bistability are 
qualitatively different for (b) on-site and (c) inter-site loca- 
tions of the resonator along waveguide and this cannot be 
distinguished within the conceptual framework of structure 
of type (a). 



ties and, as a result, may significantly enhance the ef- 
fective coupling between short pulse and resonators, and 
(ii) photonic crystals allow the creation of ultra-compact 
high-Q resonators, which are essential for the further 
miniaturization of all-optical nanophotonic devices. De- 
spite this, many researchers still believe that the basic 
properties of devices based on ridge waveguides or PhC 
waveguides are qualitatively identical, and that they can 
be correctly described by the coupled-mode theory for 
continuous systems (see Refs. [11-21] and the discussion 
in Sec. II). 

However, an inspection of Figs, l(a-c) reveals, that a 
major difference between the ridge waveguide in (a) and 
PhC waveguides in (b,c) is that a PhC waveguide is al- 
ways created by an array of coupled small-volume cav- 
ities and, therefore, exhibits an inherently discrete na- 
ture. This suggests that in these systems an additional 
coupling parameter appears which relates the position 
of the a-resonator to the waveguide cavities along the 
waveguide. As a matter of fact, we may (laterally) place 
the a-resonator at any point relative to two succesive 
waveguide cavities., thus creating a generally asymmet- 
ric device which (in the nonlinear transmission regime) 
should exhibit the properties of an optical diode, i.e., 
transmit high-intensity light in one direction only. This 
is an intriguing peculiarity of photonic-crystal based de- 
vices which we will analyze in a future publication. In 
this paper, however, we restrict our analysis to symmetric 
structures and study the cases of either on-site coupling of 
the a-resonator to the PhC waveguide, shown schemat- 
ically in Figs. 1(b), or inter-site coupling, as shown in 



Fig. 1(c). 

To address these issues, we employ a recently devel- 
oped approach [22-24] and describe the photonic-crystal 
devices via effective discrete equations that are derived by 
means of a Green's function formalism [25-29]. This ap- 
proach allows us to study the effect of the discrete nature 
of the device on its transmission properties. In particular, 
we show that the transmission depends on the location 
of the resonance frequency uj a of the a-resonator with 
respect to the edges of the waveguide passing band. If 
uj a lies deep inside the passing band, all devices shown in 
Figs, l(a-c) are qualitatively similar, and can adequately 
be described by the conventional coupled-mode theory. 
However, if the resonator's frequency uj a moves closer 
to the edge of the passing band, standard coupled-mode 
theory fails [30]. More importantly, we show that in 
this latter case the properties of the devices shown in 
Figs. 1(b) and Fig. 1(c) become qualitatively different: 
light transmission vanishes at both edges of the passing 
band, for the cases shown in Fig. 1(a) and Fig. 1(b), but 
for the case shown in Fig. 1(c) it remains perfect at one 
of the edges. Moreover, the resonance quality factor for 
the structure (c) grows indefinitely as uj a approaches this 
latter band edge, accordingly reducing the threshold in- 
tensity required for a bistable light transmission. This 
permits to achieve a very efficient all-optical switching in 
the slow-light regime. 

The paper is organized as follows. In Sec. II we sum- 
marize and extend the results of standard coupled-mode 
theory which accurately describes the system shown in 
Fig. 1(a). Then, in Sec. III. A we derive a system of ef- 
fective discrete equations [25, 26] and utilize a recently 
developed approach for its analysis [22, 23]. Specifically, 
in Sec. III.B and Sec. III.C, respectively, we study the 
two geometries of the waveguide-resonator coupled sys- 
tems schematically depicted in Fig. 1(b) and Fig. 1(c). 
In Sec. IV, we illustrate our main findings for several 
examples of optical devices based on a two-dimensional 
photonic crystal created by a square lattice of Si rods. Fi- 
nally, in Sec. V we summarize and discuss our results. For 
completeness as well as for justification of the effective 
discrete equations employed, we include in Appendix A 
an analysis of simpler cases of uncoupled cavities and 
waveguides. The effects of nonlocal waveguide dispersion 
and nonlocal waveguide-resonator couplings are briefly 
summarized in Appendix B. 



II. COUPLED-MODE THEORY 

In this Section, we first summarize the results of stan- 
dard coupled-mode theory and other similar approaches 
developed for the analysis of continuous- waveguide struc- 
tures similar to those displayed in Fig. 1(a). Then, we 
extend these results in order to obtain analytical formu- 
las for the description of bistable nonlinear transmission 
in such devices. 
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A. Linear transmission 

Transmission of light in waveguide-resonator systems 
is usually studied in the linear limit using a coupled-mode 
theory based on a Hamiltonian approach. This approach 
has been pioneered by Haus and co-workers [11, 12] and 
is similar to that used by Fano [13] and Anderson [14] for 
describing the interaction between localized resonances 
and continuum states in the context of an effect which 
is generally referred to as "Fano resonance". For the 
analysis of the transmission of photonic-crystal devices, 
this approach has been employed first by Fan et al. [15] 
and has been elaborated on by Xu et al. [16]. 

Throughout this paper we consider the propagation of 
a monochromatic wave with the frequency u) lying inside 
the waveguide passing band; we assume that the waveg- 
uide is single-moded as well as that the resonator a is 
non-degenerate and losses can be neglected. In this case, 
the complex transmission and reflection amplitudes, t(u>) 
and r(w), can be written in the form 



t(w) 



a(w) 



r{uj) 



(1) 



<r{bj) — i ' v 7 

with a certain real- valued and frequency-dependent func- 
tion a(ui) and the reflection phase <p r (w)- Accordingly, 
the absolute values of the transmission coefficient T — 
\t\ 2 and reflection coefficient R = 



are 



T(w) 



and R(u>) 



(2) 



CT 2 H + 1 " v ~ ; <7 2 H + 1 ' 

and it is easy to see that T + R = 1 for any <t(uj). 

If the frequency u a of the resonator a lies inside the 
waveguide passing band, Fano-like resonant scattering 
with zero transmission at the resonance frequency w rcs , 
lying in the vicinity of the resonator's frequency, u) a , 
should be observed [13, 22]. This corresponds to the 
condition (j(w res ) = and, based on the terminology de- 
veloped in Refs. [17, 18], o-(co) may be interpreted as the 
detuning of the incident frequency from resonance. 

The results of standard coupled-mode theory analysis 
(for instance, see Ref. [16]) indicate that in the vicinity of 
a high-quality (or high-Q) resonance, the detuning func- 
tion (j{uu) can be accurately described through the linear 
function 



- u 



7 



where 7 : 



2Q 



(3) 



which leads to a Lorcntzian spectrum. Here, Q is the 
quality factor of the resonance mode of the a-resonator. 
From the Hamiltonian approach [16], we find that the res- 
onance frequency aj ros almost coincides with the resonator 
frequency u a (see, however, Appendix A in Ref. [19] for 
a more accurate estimate of w ros ), the reflection phase is 
ip r = 7r/2, and the resonance width 7 is determined by 
the overlap of the mode profiles of waveguide and res- 
onator: 



L 



7: 



dr5e{r)£* k {r)£ a {r) 



(4) 



Here, £ a (r) is the normalized dimcnsionless electric field 
of the resonator mode, £ k {r) is the corresponding field 
of the waveguide mode at wavevector k = k(u) Tes ), v gr = 
(dio/dk) is the group velocity calculated at the resonance 
frequency, and L is the length of the waveguide section 
employed for the normalizing the modes to 

J dre wg (r)\£ k (r)\ 2 = W k , 

wg section 

Jdr s a (r)\£ a (r)\ 2 = W a . 

all space 



(5) 



Furthermore, e a (f) and e wg (f) are the dielectric func- 
tions that describe the resonator a and waveguide, re- 
spectively. From Eqs. (4)-(5) it is easy to see that the 
resonance width 7 does not depend on the length L. 

However, within the Hamiltonian approach, the func- 
tion Se(r) in Eq. (4) remains undetermined. Generally, 
it is assumed to be a difference between the total dielec- 
tric function and the dielectric function So(r) "associated 
with the unperturbed Hamiltonian" [16] which is an ill- 
defined quantity. A different approach based on a per- 
turbative solution of the wave equation for the electric 
field [20] sheds some light on the resolution of this am- 
biguity and shows explicitly that eo(r) can be taken as 
either s wg (f) or e a (r). 



B. Nonlinear transmission 

If the resonator a is made of a Kerr-nonlinear ma- 
terial, increasing the intensity of the localized mode of 
the resonator leads to a change of the refractive index 
and, accordingly, to a shift of the resonator's resonance 
frequency. As a result, the nonlinear light transmission 
in this case is described by the same Eqs. (l)-(2), with 
the only difference that the frequency detuning parame- 
ter <j{uj) should be replaced by the generalized intensity- 
dependent frequency detuning parameter (<j(uj) — J a ). 
Here, J a is a new dimensionless parameter which is, as 
we show below, proportional to the intensity of the res- 
onator's localized mode. In particular, Eqs. (2) take the 
form 



Hoj) - J a } 2 + 1 



R = 



Hui) - j a ] 2 + 1 



(6) 



In order to find an explicit expression for J a , we as- 
sume that: (i) The dimensionless mode profiles £ a (r) 
and £ k {r) introduced in Eqs. (4)-(5) are normalized to 
their maximal values (as functions in real space), i.e., 

I^WlLx = tolLx = 1; (ii) The physical electric 
fields are described by amplitudes, A a and A k , multiply- 
ing the field profiles. Consequently the maximum inten- 
sity of the electric field in the vicinity of the a-resonator, 
E(r) ~ A a £ a (r), is equal to |^4 Q | 2 ; (iii) The a-resonator 
is made of a Kerr-nonlinear material with the nonlinear 
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(3) 

susceptibility Xa and it covers the area described by the 
function 6 a {r). This function is equal to unity for all r 
inside the cavities which form the resonator structure and 
vanishes outside. In this case, J Q takes the form 



_ 12nQ K (3) 



(7) 



where k is the dimensionless and scale-invariant nonlin- 
ear feedback parameter (first introduced in similar form 
in Refs. [17, 18]) which measures the geometric nonlinear 
feedback of the system. It depends on the overlap of the 
resonator's mode profile with spatial distribution 9 a (r) 
of nonlinear material according to 



3 

W 2 



c 



dr 9 a (f)e a (r)\£ a (f)\ 4 



(8) 



all space 



where d is the system dimensionality. 

The dependence of J a on the power of the incoming 
light has already been studied analytically in Refs. [18, 
20, 21]. Here, we suggest a simpler form for this depen- 
dence 



Jin - Ja ([<T(W) - J a f + 1) 



(9) 



where we have introduced the dimensionless intensity Jj n 
which is proportional to the experimentally measured 
power of the incoming light 



Pin = 



c 2 k(m) 
2ttuj 



I'm — PaJin 



(10) 



In this expression, we have abbreviated the incoming 
light intensity as I ln = \Ak\ 2 and introduced the char- 
acteristic power Pq of the waveguide defined as (see 
Refs. [17, 18, 20, 21] for derivation): 



Po = 



d-l 



Q 2 K oXa 



(3) 



(11) 



Finally, the outgoing light power P ont = Po Jout can be 
determined through the dimensionless intensity of the 
outgoing light J out = T Ji n with the transmission co- 
efficient T defined by Eq. (6). 

It follows from Eqs. (6) and (9) that the nonlinear 
transmission problem is completely determined by the 
value of a(ui) and the sign of the product a(cu) • J Q . As is 
illustrated in Fig. 2, for frequencies where (o~(u>)- J a ) < 0, 
the transmission coefficient T and the outgoing light in- 
tensity J ou t grow monotonically with Ji n for all values 
of a(u>). 

The situation becomes more interesting for frequencies 
lying on the other side of the resonance where (cr(w) • 
Ja) > 0. In this case T and, therefore, J ou t become non- 
monotonic functions of J n , as is illustrated in Fig. 3. 
Moreover, for o~ 2 (ui) > 3 these functions become multi- 



valued functions of Ji n in the interval J 



(3,4) 



< Jin < 



J (J. ) 

fY V 111' 




FIG. 2: (Color online) Dependencies of the transmission co- 
efficient, T, the outgoing light intensity, J out , and the res- 
onator's mode intensity, J a , on the incoming light intensity, 
J7in, for <j 2 (<jj) = 1 and negative product (cr(a>) • Ja). 



J:^' 2 \ where 



J 



(1,2) 



7 



(3,4) 



27 
2 

27 



{cT 3 +9cT+[a 2 -3] 3/2 } , 

{<7 3 + 9<t- [<r 2 -3] 3/2 } , (12) 



which are also shown in Fig. 4. In this interval the nonlin- 
ear light transmission becomes bistable: low- and high- 
transmission regimes coexist at the same value of the 
incoming light intensity J mi as can be seen in Fig. 3 
for cr 2 > 3 (intermediate parts of the curves correspond 
to unstable transmission). Therefore, by increasing an 
intially low intensity Ji„ we obtain a hysteris where we 
jump from the point (1) to (2), and then upon decreasing 
Jin, we jump from the point (3) to (4). The transmission 
coefficients at these characteristic points are 



y(l,3) 
T (2,4) 



1 



2ct2(i t ^i_ 3 / (T 2)_2 
(IT 2^1-3/(7 



2\2 



5-3/ct 2 t4v/1-3/ct 2 



(13) 



and they are depicted in Fig. 4. For completeness, we also 
present the expressions for the resonator's mode intensity 
at these points 



2ct 



Ji 2 < 4) = yifv^s 



(14) 



From a practical point of view, these solutions have im- 
portant consequences. Firstly, the bistability condition 
cr 2 > 3 corresponds to a linear transmission T > 3/4. 
That is, the bistable transmission becomes possible only 
for frequencies where (<r(u)) • J a ) is positive and linear 
transmission exceeds 75%. As demonstrated in Fig. 4 
and Eq. (12), when a 1 grows, all threshold intensities 
grow, too, starting with the minimum threshold inten- 
sity J in 1,2 ' 3 ' 4) = 8/3 15 « 1.54 at a 2 = 3. 

For ideal nonlinear switching the coefficients and 
T< 4 ) should be close to unity while and should 
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FIG. 3: (Color online) Dependencies of the transmission co- 
efficient, T, the outgoing light intensity, Jout, and the res- 
onator's mode intensity, J a , on the incoming light intensity, 
J ln , for several different values of <y 2 (uj) and positive product 
{o{lo) ■ J a ). 



vanish. However, as can be seen from Fig. 4 and the 
asymptotic (for large a 2 ) expressions 

T-(i) ~ 1 - — T< 2 ) w 1 — 

' r2 ' 4(T 2 ' 



T( 3 ) « 1-1 w — 

a 2 ' 4a 2 



(15) 



of Eqs. (13), these conditions cannot be satisfied simul- 
taneously. In particular, the transmission coefficient 
does not vanish but approaches unity for large a 2 . More- 




FIG. 4: (Color online) Dependencies of the threshold incom- 
ing light intensity and the corresponding transmission coeffi- 
cients on u 2 (uj) for the four critical points (l)-(4) indicated 
by circles in Fig. 3. Here, we assume that (<j(uj) - Jo) > 0. 



over, there exists no condition under which and 
vanish simultaneously. Therefore, it is impossible to cre- 
ate ideal nonlinear switches in these systems. 

A reasonable compromise for realistic nonlinear switch- 
ing schemes of this type could be the usage of the fre- 
quency with cr 2 ~ 5, for which the linear light transmis- 
sion is close to 83%. For this case, the critical trans- 
mission coefficients T< 2 ) ~ 3.7% and ~ 7% are suf- 
ficiently small, while ~ 60% and ~ 74% are 
large enough for practical purposes. The threshold in- 



tensities J- 



1,2) 



2.53 and J, 



(3,4) 



2.11 differ about 



20% from each other, so that in this case one can achieve 
a high-contrast and robust switching for sufficiently small 
modulation of the incoming power. 

The above analysis suggests that the optimal dimen- 

fi) 

sionless threshold intensities are fixed around J>J ~ 2.5 
so that the real threshold power of the incoming light, 

^in = ^bi^in i can only be minimized by minimizing the 
characteristic power, Pq, of the system. An inspection 
of Eq. (11) shows that this can be facilitated by increas- 
ing the resonator nonlinear feedback parameter, K a , the 

(3) 

material nonlinearity, \a , or the resonator quality fac- 
tor, Q. For small-volume photonic crystal resonators, it 
has been established that k ~ 0.2 (see [17, 18]), and this 
value can hardly be further increased. 

Therefore, only two practical strategies remain that 
could lead to an enhancement of nonlinear effects in this 
system. The first approach is based on specific mate- 
rial properties: We should create the resonator a from a 

(3) 

material with the largest possible value of Xa ■ In high- 
index semiconductors, nearly instantaneous Kerr nonlin- 
earity reaches values of ri2 ~ 1.5 x 10~ 13 cm 2 /W [31], 
where ri2 oc x^ / n o an d «o is the linear refractive index 
of the material. Even such relatively weak nonlinearity 
is already sufficient for many experimental observations 
of the bistability effect in the waveguide-resonator sys- 
tems [9, 10]. However, using polymers with nearly in- 
stantaneous Kerr nonlinearity of the order of ri2 > 10~ n 
cm 2 /W and, at the same time, sufficiently weak two- 
photon absorption [32], one could potentially decrease 
the value of Pq by at least two orders of magnitude. Poly- 
mers, however, have a low refractive index which is insuf- 
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ficient for creating a (linear) photonic bandgap required 
to obtain good wavcguiding and low losses. The solution 
to this could be the embedding of such highly nonlin- 
ear but low-index materials into a host photonic crystal 
made of a high-index semiconductor. Optimized waveg- 
uding designs for the basic functional devices of this kind 
are available [33-35] and recent experimental progress 
[36, 37] may soon allow a realization of corresponding 
linear and nonlinear devices. 

The second approach is based on designing waveguide- 
resonator structures with the largest possible quality fac- 
tor, Q. Potentially, one can increase Q indefinitely by 
mere increase of the distance between the waveguide and 
the resonator. However, this leads to a corresponding 
increase in the size of the nonlinear photonic devices. A 
very attractive alternative possibility for increasing Q is 
based on the adjustment of the resonator geometry [38]. 

In what follows, we suggest yet another possibility to 
dramatically increase Q through an optimal choice of the 
resonator location relative to the discrete locations of the 
cavities that form the photonic-crystal waveguide. 



C. Limitations of the coupled- mode theory 

Standard coupled-mode theory exhibits a number of 
limitations. Firstly, it gives analytical expression for the 
detuning parameter a(uj) only near the resonator fre- 
quency uj a . And this immediately highlights the second 
limitation: standard coupled-mode theory [16-21] can- 
not analytically describe resonant effects near waveguide 
band edges. However, numerical studies [30] have re- 
cently demonstrated that the effects of the waveguide 
dispersion become very important at the band edges and 
may lead to non-Lorentzian transmission spectra in cou- 
pled waveguide-resonator systems. 

As a matter of fact, the question "what happens if the 
resonator frequency cj a lies near the edge of the waveg- 
uide passing band or even outside it?" may be of a great 
practical importance due to two reasons. Firstly, in real- 
istic structures it is not always possible to appropriately 
tune the frequency u> a , and therefore it is important to 
understand properties of the system for any location of 
the resonance frequency. Secondly, as we have already 
mentioned in the Introduction, PhC waveguides can pro- 
vide us with a very slow group velocity of the propagating 
pulses — but in most cases they do it exactly at the pass- 
ing band edges. Therefore, if we wish to utilize such a 
slow light propagation for nonlinearity enhancement, we 
should extend the above analysis to such cases, too. 

In what follows, we describe an alternative analyti- 
cal approach to the coupled waveguide-resonator struc- 
tures which allows us to correctly analyze both linear 
and nonlinear transmission for arbitrary locations of the 
resonator frequency uj a relative to the waveguide passing 
band, including the transmission near band edges in the 
slow light regime. 



III. DISCRETE MODEL APPROACH 

Having discussed the results obtained for the 
continuous- waveguide structure shown in Fig. 1(a), we 
now take into account the discrete nature of the waveg- 
uding structure embedded in photonic crystals. In par- 
ticular, we analyze what will change in the system prop- 
erties when we move the resonator along the waveguide 
from the on-site location shown in Fig. 1 (b) to the inter- 
site location shown in Fig. 1(c). Our analysis is based on 
effective discrete equations that have been derived for the 
description of photonic crystal devices [25-29] in combi- 
nation with a recently developed discrete model approach 
to nonlinear Fano resonances [22]. 



A. Discrete Equation Approach 

First, we derive an appropriate set of discrete equations 
[see Eqs. (24) below], and show that they can be applied 
to a variety of the photonic-crystal devices. We start 
from the wave equation in the frequency domain for the 
electric field 



V x V x 



o 2 ^ 



E{r) = 



(16) 



where the dielectric function e(r) = i pc (r) + Se(f) con- 
sists of the dielectric function i pc (r) of a perfectly peri- 
odic structure and a perturbation 5e(r) that describes the 
embedded cavities. It is convenient to introduce the ten- 
sorial Green function of the perfectly periodic photonic 
crystal, 



VxVx-(-) e pc (r) 



G(f,r» = IS(f-r') (17) 



and to rewrite Eq. (16) in the integral form, 

E(r) ={^Y J r' G(r,r'\u)Se{r')E(r') , (18) 

where we assume that the frequency ui lies inside a com- 
plete photonic bandgap so that the electric field vanishes 
everywhere except for areas inside and in the vicinity of 
cavities. We enumerate the cavities by an integer index 
n and introduce dimensionlcss functions 9 n (f) which de- 
scribe the shape of the n-th cavity. As a result, Se(r) 
may be represented as 



Se(r) 



E 



5e n + xi 3) |-B(r)| 2 n (f-R n ), (19) 



where R n , Se n , and Xn^ are, respectively, position, (lin- 
ear) dielectric function, and nonlinear third-order suscep- 
tibility of the n-th cavity. 

Similar to Sec. II, we describe the electric field of the 
n-th cavity mode via a dimensionless field profile £ n (r) 
and a complex amplitude A n . Taking into account that 
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inside the cavities the electric field of the system is a 
superposition 



(20) 



Eq. (18) can be rewritten as a set of discrete nonlinear 
equations 



D„{u)A r , 



V n , m (w)A n 



«„Hxl 3) |A n | 2 A„,(21) 



where D n {uj) = 1— V n , n (u>) is the dimensionless frequency 
detuning from the resonance frequency, ui n , of the n-th 
cavity. Furthermore, 

V n . m {uj) = (~Y / 'dr I dr' £*{f)e n {f) (22) 



W n Vc- 

x 6 m (r')G(r + R n - R m , r » £ m (r ') , 

is the dimensionless linear coupling between the n-th and 
the m-th cavity. Similarly, 

Kn{Uj) = W n OY ^ (23) 

x e n (r')G(r,r'\Lu)\£ n (r')\ 2 £ n (r') , 

is the dimensionless and scale-invariant nonlinear feed- 
back parameter which should be compared with the 
analogous parameter (8) introduced in the conventional 
coupled-mode theory analysis [17, 18]. Finally, W n is 
defined in exactly the same way as W a in Eq. (5). 

We remark that in deriving Eqs. (21) we have neglected 
higher-order couplings proportional to the integrals of 
£*(f)£ m (r + R n — R m ) with n ^ m but take into ac- 
count the coupling coefficients which involve integrals of 
G(f + R n — Rm,r'\uj) with n ^ m. This approxima- 
tion is sufficiently accurate in most cases, as we demon- 
strate in Refs. [24, 26]. We would like to mention that 
in Eqs. (21)-(23) we have used more accurate definitions 
of the coupling coefficients than those that have been in- 
troduced earlier in Refs. [26-29]. They have also a more 
generic form than those we used in Ref. [25]. 

Typical frequency dependencies of the parameters of 
the discrete model, Eq. (21), are displayed in Figs. 9- 
11 of Appendix A, where we also discuss the application 
of Eqs. (21)-(23) to simple structures such as linear and 
nonlinear photonic crystal resonators and straight waveg- 
uides. Here, we apply Eqs. (21)-(23) to study the more 
complicated case of the nonlinear coupled waveguide- 
resonator systems shown in Figs. l(b,c). The set of 
Eqs. (21) may be separated in this case according to 

L 

D w (w)A n = ^ V.j W (uj)(A n+j + A n _j) + V nta (w)A a , 

(24) 



where we assume that all cavities of the photonic-crystal 
waveguide are identical and linear, so that we can de- 
note D w (lu) = D n (uj) and Vj w (uj) = V ntn ±j{w) for any n 
inside the waveguide. Furthermore, the index a defines 
the parameters of the side-coupled nonlinear resonator. 
Below we show that the assumption of linear waveguide 
cavities may be relaxed for frequencies near the resonator 
resonance frequency iv a because then the amplitudes A n 
remain small in comparison with the amplitude A a . 

For the first equation in Eq. (24), we seek solutions of 
standard form 



r l/2 



An. — 



ik(uj)sn 



for n 3> 1, 



jl/2 el k(u)sn + r ^ e -tk(u)sn for n < _^ 



,(25) 



where s is the distance between the nearest waveguide 
cavities and I{ n is the intensity of the incoming light. 
For both structures shown in Figs. l(b,c), we obtain that 
the transmission and reflection coefficients can formally 
be described by the same expressions (l)-(2) as for the 
structure depicted in Fig. 1(a). However, within the dis- 
crete equation approach the expression for the detuning 
parameter o~{u) can now be found for the entire frequency 
range. Below, we discuss novel results for the structures 
shown in Fig. 1(b) and Fig. 1(c) separately. 



B. On-site resonator 

First, we obtain the solution of this problem for the 
structure shown in Fig. 1(b). For simplicity, we as- 
sume that the only nonvanishing coupling coefficients in 
Eq. (24) are Vi w (cu), V a fi(w), and V~o. a (uj) (see, however, 
Appendix B for a more accurate analysis which takes into 
account additional coupling coefficients). As a result, 
we obtain the transmission and reflection coefficients de- 
scribed by Eqs. (l)-(2) with <p r = ir/2 and a correspond- 
ing expression for a(u)): 



<t(lj) = 2 sin [fc(w)s] 



ViwM A 

V , a (w) A a 



(26) 



which should be considered as a generalized intensity- 
dependent frequency detuning parameter <r(ui) + J a in- 
troduced in Eq. (6) above. The amplitude Aq in Eq. (26) 
is given by 



A = t(w)L, 



1/2 



(27) 



while the waveguide dispersion relation k{u>) is deter- 
mined by Eq. (A6). 

In the case of a linear a-resonator (i.e. = 0), 

the amplitude A a = V a fi(<-<j)Ao/ D a (u>) is proportional to 
the amplitude A . Therefore, <j(u>) and, accordingly, the 
transmission and reflection coefficients do not depend on 
the light intensity. Upon introducing the abbreviation 



D a (uj)Vi w (uj) 
Vo, a (u)V at o(v) 



(28) 
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0.38 0.39 0.4 
Frequency, a/A, = (oa/27tc 



FIG. 5: (Color online) Linear transmission through a photonic 
crystal waveguide that is created by removing every second 
rod in a row (s = 2ai) side-coupled to a one-site resonator cre- 
ated by removing a single rod. The underlying 2D photonic 
crystal is described in Appendix Al. We compare exact nu- 
merical results (solid line) with the analytical results based 
on Eq. (29) (dot-dashed line) and Eq. (Bl) (dashed line). 



the detuning parameter, Eq. (26), for a linear a-resonator 
reads as 



cf(lu) — 2fi(u>) sin [fc(w)s] 



(29) 



This implies that cr(w) vanishes when either D a (uj) = 
or k{uj) = Ttn/s with an arbitrary integer n. The first 
condition reproduces Eq. (3) with uj tcs = uj a and the 
resonance width 7 given by 



S0J a Ul v 



7 



' V a V w Vo,q V a .o 



sm[fc(w a )s] 1 
where v a and v w are defined by Eq. (Al), 

Vo. a V a n Vo. a V a n 



2u) a D' a Vx 
and the group velocity 

dk 



2V lv 



2s w w i/ w Tq w sin[k(L0 a )i 



(30) 



(31) 



(32) 



can be found directly from Eq. (A6). Here and in what 
follows, we assume that the values of all frequency- 
dependent parameters whithout explicitly stated fre- 
quency dependence are evaluated at the resonance fre- 
quency, lu ics . Finally, we notice that the resonance width, 
Eq. (30), is very similar to that described by the coupled- 
mode theory, Eq. (4). 

It is important that the quality factor Q of the reso- 
nance 



Q = 



uj a sin [fc(aj Q )s] v g 



1 



27 



2A„ 



suj w 2 v a v w Vq a V a 



, (33) 



is multiplied by the factor sin [k(uj a )s] ~ v gr , and, there- 
fore, becomes strongly suppressed near the edges of 
waveguide passing band, k(ui a ) = 0, ±ir/s. Accordingly, 
the detuning parameter (29) vanishes at these edges, too. 
This means that, in agreement with the numerical calcu- 
lations shown in Fig. 5, the transmission coefficient T{w) 
vanishes not only at the resonance frequency, but also at 
both edges of the waveguide passing band. Such an effect 
was recently observed by Waks and Vukovic [30] in their 
numerical calculations based on standard coupled-mode 
theory which takes into account the waveguide disper- 
sion. Therefore, the effect of vanishing transmission at 
the spectral band edges may be attributed also to the 
structure shown in Fig. 1(a). 

Obviously, this enhancement of light scattering at the 
waveguide band edges should be very important from the 
point of view of fabrication tolerances since virtually any 
imperfection contributes to scattering losses. Moreover, 
as discussed in Sec. IV, this effect is detrimental to the 
concept of all-optical switching devices based on slow- 
light photonic crystal waveguides. 

We support this conclusion by another observation. 
First, the light intensity at the 0-th cavity, \Aq\ 2 = 
T(w)ij n , vanishes at the resonance frequency for arbitrary 
large incoming light intensity, because T(uj a ) = 0. There- 
fore, the nonlinearity of this cavity may safely be ne- 
glected. In contrast, the light intensity at the a-resonator 
reaches its maximum value at uj a , 



A a (u a )\ 

He 



2^4' V * 



V , a 



sin 2 [k(cj a ) s] ■ /j, 



I ia ^(2Qu a V a , y-I- m , (34) 



which may significantly exceed the incoming light inten- 
sity 7i n when the coupling Vo. a between the a-resonator 
and waveguide becomes small enough relative to the cou- 
pling Vi w between the cavities in the waveguide. This 
strong enhancement suggests a physical explanation for 
the existence of the rather strong nonlinear effect of light 
bistability at relatively low intensities of the incoming 
light. However, when the resonance frequency uj a lies 
close to any of the waveguide band edges, it is seen from 
Eq. (34) that the light intensity at the a-resonator be- 
comes (strongly) suppressed by a factor sin 2 [fc(a; Q )s] . 

Details of an extension of the above discussion to the 
case of more realistic non-local couplings, i.e., more than 
nearest neighbors couplings, is presented in Appendix B 
and here we only summarize the results. Both, a non- 
locality of the inter-coupling between waveguide cavities 
as well as a nonlocality of cross-coupling with the a- 
resonator lead to a small shift in the resonance frequency, 
u> Ies , but do not change the main result about the sup- 
pression of the detuning a(co) and transmission T(u>) at 
both edges of the waveguide passing band. However, we 
would like to emphasize that for a fully quantitative anal- 
ysis, non-local couplings have to be taken into account, 
for instance, within the framework of the recently devel- 
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oped Wannier function approach [39]. 

We now consider the case when the resonator a is non- 
linear, i.e. \a ^ 0. As has been previously shown 
in Ref. [22], this case, too, can be studied analytically 
even for non-local couplings between the cavities and 
resonator and novel effects originating solely from the 
non-locality may be expected when the non-local cou- 
pling strength exceeds one half of the local coupling. Un- 
fortunately, in realistic photonic crystals this limit may 
hardly be realized so that here we restrict our analysis to 
the local-coupling approximation. In this case, we obtain 
from the second equation in Eqs. (24) that the amplitude 
A a uniquely determines the amplitude Aq. Substituting 
the latter expression into Eqs. (26)-(27), we find that 
the nonlinear transmission is described by Eqs. (6) and 
(9) with the detuning cr(u)) determined by Eqs. (28)-(29) 
and the dimensionless intensities J a and J- m given by the 
expressions 



(35) 



6e 



8sin[fc(w rcs )s] V lw Q 2 K Q ^ x (3)/ in 

slu w \de a v w J 

where Q is determined by Eq. (33). Therefore, all the 
results for the nonlinear light transmission which are dis- 
played in Figs. 2-4 are directly applicable to the structure 
of Fig. 1(b), too. 

In an experiment, one measures not the light inten- 
sity in the waveguide, ii n , but the propagation power, 
Eq. (10), where for the discrete structure of Fig. 1(b), 
the characteristic power Pq is 



Po 



C 2 k(u> a ) 



167r sm[k(u> a )s\u} a Vi w 

C 2 k(uj a )s f LJwfec^w 



Se 



1 



8nv e 



UJaSeQVa 



1 



(3) 



(3) 



(36) 



Again, this result is quite similar to Eq. (11) for the con- 
tinuous structure of Fig. 1(a). Nevertheless, our more 
general analysis explicitly suggests that it should be bet- 
ter to use the a-resonator with the resonance frequency at 
the center of the waveguide passing band k{uj a ) w tt/2s, 
where the group velocity reaches its maximum. No- 
tice, however, that this suggestion becomes wrong for 
the structure of Fig. 1(c) studied in the next subsection. 



C. Inter-site resonator 

In the system where the a-resonator is placed sym- 
metrically between two cavities of the waveguide and, 
therefore, couples equally to both of them, a qualitatively 
different type of resonant transmission occurs. The cor- 
responding structure is schematically shown in Fig. 1(c). 
Assuming that in this case the nonvanishing coupling co- 
efficients in Eq. (24) are Vi w (u), V a ,i(u)) = V at o(w), and 




0.38 0.39 0.4 0.41 
Frequency, a/A, = ooa/27tc 



FIG. 6: (Color online) Liner transmission through a photonic 
crystal waveguide created by removing every second rod in a 
row (s = 2ai) side-coupled to an inter-site resonator created 
by removing a single rod. The underlying 2D photonic crystal 
described in Appendix A 1. We compare exact numerical re- 
sults (solid line) with the analytical results based on Eq. (39) 
(dashed line). 



Vi t0 ,(uj) = Vo,a(w), we seek solutions to the first equa- 
tion of the system (24) that are of the form of Eq. (25). 
Again, we find that the transmission and reflection coeffi- 
cients are given by Eqs. (l)-(2) albeit with the frequency- 
dependent phase 4> r {u) = 7r/2 + k(uj)s. Here, k(uj) is 
determined by Eq. (A6), and the generalized intensity- 
dependent frequency detuning is 

7-1/2 

a{uj) + J a = i-i(e lk ^ s -l 



Mi. 

Vo, a (w) A a 



(37) 



The corresponding amplitudes are 



A - T 1 ' 2 
A - An - 



1 



[1- 

ik(ui)s jl/2 



A a 



ife M s ] Fiw(w) 

1 V Q , a (u) 



A a . (38) 



[l-e-*M ! ] V 1w (uj) 

Despite the complex form of Eq. (37), we would like to 
emphasize that the detuning a{u) determined by Eq. (37) 
is a real-valued function (see also the discussion following 
Eq. (1) above). 

In the case of the linear a-resonator (i.e., for \ a = 0), 
we obtain 

' k((jj)s s 



a(oj) = [1 + h(uj)] tan 



(39) 



where fJ-(uj) is given by Eq. (28). For a high-quality a- 
resonator in the vicinity of the resonance frequency this 
detuning parameter can be approximated by Eq. (3) with 
lu ics ~ uj a (1 — 2A Q ) and 

= 

tan[fc(w rcs )s/2] 

Here, A a is defined by Eq. (31). In contrast to Eq. (33), 
the corresponding quality factor 

tan [fc(u ros )s/2] V\ w tan [fc(w res )s/2] 



Q = 



ul rc 



27 



4A r 



2 v a Vo.qV^o 



L(41) 
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is now multiplied by the factor tanJfes/2] which does not 
vanish and even diverges as k(uj ros ) approaches the edge 
of the transmission band k — ±ir/s. At this band edge, 
<j(lu) ~ tan(fc(w)s/2) — > oo and, therefore, light trans- 
mission is always perfect. This conclusion is supported 
by the exact numerical calculations presented in Fig. 6. 
At the other band edge, i.e., for k(uj) = 0, transmission 
vanishes, similar to the structures shown in Figs. l(a,b). 

The light intensity at the a-resonator reaches its max- 
imal value at the resonance frequency 



\A a (tO Tes )\' 



sin 



V , a 

4:QlS a V a} COS 



fc(cj rcs )s 
2 

k(uj TCS )s 



■In 

2 



Iin (42) 



Again, in contrast to the corresponding light intensity 
(34) for the on-site coupled structure, Eq. (42) does not 
vanish at the edge of the transmission band k = ±ir/s. 
Therefore, we can expect that for inter-site coupled struc- 
ture nonlinear effects at the band edge k = ±tt/s should 
be sufficiently strong to allow bistable transmission and 
switching. 

To investigate this, we assume that the a-resonator 
is nonlinear (x a 0) an d introduce the same dimen- 
sionless intensities J a and J ln as in Eq. (35). However, 
now the quality factor Q is defined by Eq. (41) and the 
resonance frequency is to ICS ~ uj a (l — 2A Q ). We find 
that this nonlinear problem, too, has a solution of the 
form given by Eqs. (6) and (9). However, now the de- 
tuning u{uj) is given by Eq. (39). Therefore, all results 
presented above in Figs. 2-4 remain applicable to this 
structure, too. The only but very important qualitative 
difference of the structure shown in Fig. 1(c) is that the 
transmission coefficient T(uj) and the corresponding light 
intensity \A a \ 2 at the a-resonator do not vanish at the 
band edge k — ±7r/s since the quality factor Q at this 
band edge grows to infinity for the inter-site structure of 
Fig. 1(c). Therefore, this structure may be utilized for 
realizing efficient all-optical switching devices based on 
slow-light photonic crystal waveguides. This is in sharp 
contrast to the structures shown in Figs. l(a,b). 



IV. DISCUSSION OF RESULTS 

In this section, we summarize our results and em- 
phasize their importance by applying them to spe- 
cific photonic-crystal structures. We consider a two- 
dimensional photonic crystal created by a square lattice 
of dielectric rods in air. The rods are made from Si or 
GaAs (e — 11.56) and have radius r = 0.18a. 

First, we consider a waveguide created by removing 
every second rod (s = 2a) in a straight line of rods cou- 
pled to a nonlinear resonator a created by replacing a 
single rod of the two-dimensional lattice with a highly- 
nonlinear polymer rod. The corresponding structure is 
schematically shown in the insets in Fig. 7. The resonant 
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FIG. 7: (Color online) Linear transmission spectrum for a 
photonic crystal waveguide created by removing every second 
rod in a row (s = 2ai) side-coupled to a single on-site (a) 
or inter-site (b) polymer-rod resonator (marked by the open 
circle in the insets). The underlying 2D photonic crystal is de- 
scribed in Appendix A 1 and results for three different values 
of the resonator dielectric constant e a axe shown. 



frequency of the polymer-rod resonator lies very close to 
the edge k = ±7r/s of the waveguide passing band, and 
can be tuned by changing the linear dielectric constant 
e a of the rod. 

In Fig. 7(a) and (b) , respectively, we display the trans- 
mission spectra for both on-site and inter-site positions 
of the side-coupled resonator for three different values of 
resonator dielectric constant e a . We notice that in the 
case of the on-site position of the resonator the trans- 
mission coefficient T(lo) remains below the critical value 
of T = 75% required for bistable switching operation 
for all frequencies uj below the resonance frequency a7 rcs . 
The condition to < uj rcs corresponds to the condition 
(cr(cj) • J a ) > which should be satisfied to realize non- 
monotonic dependencies of the nonlinear transmission 
shown in Fig. 4). Therefore, this on-site system cannot 
exhibit bistability. 

On the other hand, bistability may be realized for the 
inter-site position of the side-coupled resonator for which, 
in a full agreement with our analysis presented above, the 
transmission remains perfect at the band edge k = ±ir/s 
and the quality factor Q increases as the resonant fre- 
quency approaches this band edge. In Fig. 8(b) (exam- 
ple A) we show that in this case bistable transmission 
indeed occurs for the frequency marked by a filled circle 
in Fig. 8(a). This corresponds to T(lu) = 80%, i.e., the 
choice o- 2 (ui) — 4 for the detuning parameter. 

We want to emphasize that the large value of the qual- 
ity factor (41) for the inter-site structure at k(ui) close to 
±7r/s leads to very low bistability thresholds as compared 
to the cases of on-site coupled and continous waveguide 
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FIG. 8: (Color online) (a) Linear transmission spectrum 
and (b) nonlinear bistable transmission for three different 
side-coupled waveguide-resonator photonic crystal structures 
whose topology is shown in (c). The rods consist of Si or GaAs 
(full black circles) or polymer (open red circles). Example A 
represents a close to optimal structure with inter-site loca- 
tion of the a-resonator whose resonance frequency lies close 
to the edge k = ±7r/s of the passing band; example B rep- 
resents a sub-optimal structure with an inter-site location of 
the a-resonator whose resonance frequency lies near the cen- 
ter of the passing band; example C represents a sub-optimal 
but commonly used structure with an on-site location of the 
a-resonator whose resonance frequency lies near the center of 
the passing band. Closed circles in (a) indicate frequencies 
with T(lo) — 80% that are used for achieving high-contrast 
bistability in (b). Red circles in (c) indicate positions of the 
nonlinear polymer rods with e a — 2.56. Other parameters of 
the 2D photonic crystal are described in Appendix A 1. 



coupled structures. This is illustrated in examples B 
and C of Fig. 8: Relative to the waveguide design in ex- 
ample A, the design in example B moves the resonance 
frequency deeper into the passing band thus decreasing 
the quality factor (41). Nevertheless, the inter-site cou- 
pled example B still exhibits a much smaller bistability 
threshold than the on-site coupled system with the same 
waveguide design in example C. This is caused by (usu- 
ally) much smaller waveguide-resonator coupling and, ac- 
cordingly, much larger Q in the inter-site structures as 
compared to the on-site structures. 

Summarizing, the inter-site structure of the resonant 
waveguide-resonator interaction schematically shown in 



Fig. 1(c) allows to achieve much higher values for the 
linear quality factor Q. As a consequence, much smaller 
bistability threshold intensities for the nonlinear trans- 
mission are obtained. To employ these advantages, the 
wavevector fc(w ros ) of the guided mode at the resonance 
frequency uj ics , Eq. (39), should be as close as possible to 
7r/s. This requirement coincides with the condition of a 
very small group velocity in the waveguide and, in con- 
trast to the continuous-waveguide and on-site structures 
depicted in Figs. l(a,b), provides us with a possibility to 
create low-threshold all-optical switching devices based 
on slow-light photonic crystal waveguides. 



V. CONCLUSIONS 

We have presented a detailed analysis of PhC waveg- 
uides side-coupled to Kerr nonlinear resonators which 
may serve as a basic element of active photonic-crystal 
circuitry. First, we have extended the familiar approach 
based on standard coupled-mode theory and derived ex- 
plicit analytical expressions for the bistability thresholds 
and transmission coefficients related to light switching in 
such structures. Our results reveal that, from the point 
of view of bistability contrast (a small difference between 
two threshold intensities and robustness of switching) the 
best conditions for bistability are realized for those pa- 
rameter values for which the dimensionless detuning pa- 
rameter a(u>) is close to Vb. Practically, this corresponds 
to the choice of operation frequencies for which the linear 
light transmission is close to 83%. 

We have pointed out that the conventional coupled- 
mode theory does not allow to describe the light trans- 
mission near the band edges, and we have developed an 
improved semi-analytical approach based on the effective 
discrete equations derived in the framework of a consis- 
tent Green's function formalism. This approach is ideally 
suited for a qualitative and semi-quantitative description 
of photonic-crystal devices that involve a discrete set of 
small-volume cavities. We have shown that this novel 
approach allows to adequately describe light transmis- 
sion in the waveguide-resonator structures near the band 
edges. Specifically, we have demonstrated that while the 
transmission coefficient vanishes at both spectral edges 
for the on-site coupled structure (see Fig. 1(b)), light 
transmission remains perfect at one band edge for the 
inter-site coupled structure (see Fig. 1(c)). These fea- 
tures allow a significant enhancement of the resonator 
quality factor and, accordingly, a substantial reduction 
of the bistability threshold. As a consequence, we refer 
to this type of nonlinearity enhancement as a geomet- 
ric enhancement. The possibility of such enhancement is 
a direct consequence of the discreteness of the photonic 
crystal waveguide and is in a sharp contrast to similar 
resonant systems based on ridge waveguides. The poten- 
tial of this novel type of the nonlinearity enhancement 
may be regarded as an additional argument to support 
the application of photonic-crystal devices in integrated 
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photonic circuits. 

In addition, we would like to emphasize that the engi- 
neering of the geometry of photonic-crystal based devices 
such as that presented in Fig. 1(c) becomes extremely 
useful for developing novel concepts of all-optical switch- 
ing in the slow-light regime of PhC waveguides which may 
have much wider applications in nanophotonics and is 
currently under active experimental research [40] . 

We believe that the basic concept of the geometric en- 
hancement of nonlinear effects based on the discrete na- 
ture of photonic-crystal waveguides will be useful in the 
study of more complicated devices and circuits and, in 
particular, for various slow-light applications. For in- 
stance, this concept may be applied to the transmission 
of a side-coupled resonator placed between two partially 
reflecting elements embedded into the photonic-crystal 
waveguide where sharp and asymmetric line shapes have 
been predicted with associated variations of the transmis- 
sion from 0% to 100% over narrow frequency ranges [41]. 
Similarly, the concept can be extended to a system of 
cascaded cavities [42] and three-port channel-drop fil- 
ters [43] , optical delay lines [44] , systems of two nonlinear 
resonators with a very low bistability threshold [45], etc. 
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APPENDIX A: CALCULATION OF THE MODEL 
PARAMETERS AND EXAMPLES 

1. Coupling coefficients for two-dimensional 
photonic crystals 

To obtain deeper insight into the basic properties of 
the effective discrete equations (21), we should know how 
the coupling coefficients D n (ui), K n (w), and V n<m (u)) de- 
pend on frequency u>. As an illustration, we consider a 
two-dimensional model of a photonic crystal consisting 
of a square lattice (lattice spacing a) of infinitely long 
dielectric rods (see Refs. [15, 17, 18] and also references 
[7-16] in Ref. [25]). We study light propagation in the 
plane of periodicity, assuming that the rods have a radius 
r = 0.18a and a dielectric constant of £ ro d = 11.56 (GaAs 
or Si at the telecommunication wavelength A ~ 1.55 /im). 
For light with the electric field polarized along the rods 
(-©-polarized light), this photonic crystal exhibits a large 




FIG. 9: (Color online) Frquency dependence of the detun- 
ing coefficient D n (u)) for the 2D photonic crystal described 
in Appendix Al, for two types of resonators: Removing a 
single rod (e„ = 1.0; solid line) leads to a localized mode at 
u>„ — 0.392, while replacing a single rod by a geometrically 
identical rod made of polymer (e„ = 2.56, dashed line) leads 
to a localized mode at ui n — 0.374. 



(38% of the center frequency) photonic bandgap that ex- 
tends from w = 0.303 (27rc/a) to v = 0.444 (27rc/a) 

Our task is to evaluate the coupling coefficients D n (oj) , 
k„(w), and V n ^ m {ijj) using Eqs. (22)-(23) with the 
Green's function G(f, r'\u>) calculated by the method de- 
scribed earlier in Refs. [28, 29]. The results of these cal- 
culations are displayed in Figs. 9-11. 



2. Isolated optical resonators 

For the case of an isolated (V n ^ m = 0) linear (xn = 0) 
optical resonator at the site n, Eq. (21) takes a sim- 
plest possible form, D n (uj)A n = 0. In this case, we 
only need to know the dimensionless frequency detun- 
ing coefficient, D n (ui). In Fig. 9 we plot D n (u>) for 
two types of resonators: a resonator created by re- 
moving a single rod and a resonator created by replac- 
ing a single rod with a polymer rod of the same ra- 
dius and e n — 2.56. Introducing the dimensionless fre- 
quency u) = a/X = (ua/2'Kc), we can express these 
coefficients, with a very good accuracy in the range 
0.36 < uj < 0.41, by the following cubic dependencies: 
D n {uj) = 9.426(w-w n )-10.889(a;--w n ) 2 +840.36(w-<I; n ) 3 
with o) n = 0.3919, for the removed rod, and D n (u>) — 
9.047(£ - w n ) - 49.555(J> - Cj n ) 2 + 770.14(£ - £ n ) 3 with 
u) n = 0.3744, for the replaced rod. 

The resonator mode can only be excited at the res- 
onator frequency oj n , which is determined by the equation 
D„(w„) = 0. Fig. 9 suggests that changing the dielectric 
constant of the resonator e„ allows to tune the frequency 
Lj n . In all cases, in the vicinity of the resonator frequency 
u> n , the coupling coefficient D n (u>) can be approximately 
expanded into the Taylor series with a linear dependence 

. , UJ — UJ n 1 , . „ 

Ai(w) ~ -, v n = , Al 

v n u n w n D' n {uj n ) 



13 



3 




-0,2 



0,32 



0,34 0,36 0,38 0,4 0,42 

Frequency, a/A, = G)a/27tc 



0.44 



FIG. 10: Frequency dependence of the nonlinear feedback 
parameter K n (ui) for the 2D photonic crystal described in 
Appendix Al. The nonlinear resonator is created by re- 
placing a single rod with a polymer rod of the same radius 
which supports at e„ = 2.56 a localized mode with frequency 
Lo n = 0.374. 



0.3 



0.2 



0,1 4. 



v „=v, 

cc,0 lw 



1— ■ I i>i-U, 



0,32 0,34 0,36 0,38 0,4 0,42 0,44 

Frequency, a/X = 0)a/27ic 

FIG. 11: (Color online) Frequency dependence of the coupling 
coefficients Vn,m(w) for the 2D photonic crystal described in 
Appendix A 1 with a on-site side-coupled waveguide-resonator 
system shown in Fig. 1-b. The notations are the same as those 
in Eq. (24) and we assume that the waveguide is created by 
removing every second rod in a row (located at R n = 2Sin) 
whereas the on-site resonator is created by replacing a single 
rod at Ru = — 2S2 with a polymer rod (e n = 2.56) of the 
same radius. The dispersion relation for such a waveguide is 
displayed in Fig. 12. 



where we have introduced a dimensionless parameter v n 
which describes the resonator sensitivity to a change of 
the dielectric constant. For our example of a polymer-rod 
resonator, we find v n w 0.295. 



(3) 



When the n-th resonator is nonlinear (i.e., xti 
0), Eq. (21) reduces to the equation D n (ui)A n = 

K n (uj) \A n \ 2 A n with the new important coefficient — 
nonlinear feedback parameter K n (u>). In Fig. 10, we de- 
pict the frequency dependence of K n (oj) for the case of 
a nonlinear polymer resonator. In the frequency range 
0.36 < to < 0.41, this behavior can be approximated 
as K n (w) = -0.111 + 1.005(a) - w n ) - 5.501(6 - w n ) 2 + 
85.57(£-a>„) 3 with u n = 0.3744. Therefore, in the vicin- 
ity of the resonator's frequency, uj n , we may assume that 



K n (oj) ~ —0.111 is constant and can rewrite Eq. (21) 
according to 



\Ar. 



(uj)Xn K n (ul n )x 



(3) 



(w-w„). (A2) 



The solution of the above equation gives us the depen- 
dence of the resonator frequency uj rcs on the resonator's 
mode intensity \A n \ 2 as 



0J n (l + KnVnX^ |A„| 2 ) 



(A3) 



Here, we have used the notation n n = K n (u> n ). As we see, 
the nonlinear sensitivity of the resonator at the site n is 
a product of its nonlinear feedback parameter, K n , the 
sensitivity to a change of the dielectric constant, v n , and 

(3) 

the Kerr susceptibility of material, Xn ■ The sign of this 
product defines the direction of the resonator frequency 
shift. In particular, for the polymer resonator used in 
Figs. 9-10, we obtain a rather small shift, n n v n w —0.033 

(3) 

which indicates that for x„ > the resonator frequency 
decreases as the light intensity grows. Designing optical 
resonators with larger n n or v n , may allow to enhance 
their nonlinear properties for a given material with Kerr 



nonlinearity x 



(3) 



3. Straight waveguides 



Now let us consider an array of identical coupled 
cavities separated by the distance s = \s\ which cre- 
ate a straight photonic-crystal waveguide depicted in 
Figs. l(b,c). Before proceeding, we would like to em- 
phasize that our analysis can equally well be applied to 
the coupled-resonator optical waveguides (CROWs) sug- 
gested in Ref. [46]. If we neglect nonlinear effects (assum- 

(3) 

ing that either the waveguide cavities are linear, Xn = 0, 
or the light intensity in the waveguide remains sufficiently 
small), Eq. (21) reduces to 



D w {u))A Tl 



3 = 1 



[co)(A 



(A4) 



Here we have defined, similar to Eq. (24), D w (oj) = 
D n {oj) and Vj W (u>) = V nt n±j{u)) which are identical for 
all n. 

In Fig. 11 we plot the frequency dependencies of 
Vi w (lu) and V2 W (w) for a photonic-crystal waveguide cre- 
ated by removing every second rod in a row, either with 
s = 2d± or with s = 2di- In the vicinity of the polymer- 
rod resonator frequency, the coupling coefficients are to 
lowest order constant: Vi w « 0.096 and Viw ~ 0.0086 . 
In the general case, our calculations show that the co- 
efficients Vj W (u)) decay nearly exponentially with j. In 
terms of frequency, they take on a constant value at the 
central passing band frequency and grow rapidly towards 
the low-frequency bandgap edge. 
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Wave vector (ks/2n) 
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FIG. 12: (Color online) Dispersion relation for a photonic 
crystal waveguide created by removing every second rod in 
a row (s — 2ai) in the 2D photonic crystal described in Ap- 
pendix A 1. Numerical exact results (solid line) are calculated 
with the super-cell plane-wave method [47] and the approxi- 
mate results are obtained from Eq. (A5) with L — 1 (dotted 
line) and L — 2 (dashed line), using the coupling coefficients 
from Fig. 11. 



According to the Floquet-Bloch theorem, Eq. (A4) has 
a solution A n = Aq cxp[±ifc(cj)sn] with an arbitrary com- 
plex amplitude Aq. The corresponding dispersion k[u>) 
is determined by the equation 



D w (oj) = 2j2Vj W (^)cos [k(uj)sj] 



(A5) 



where we assume that the coupling coefficients Vj W (u>) 
vanish for all j above L. As a matter of fact, our studies 
indicate that sufficiently accurate results can be obtained 
already for L ~ 4a/ a. In Fig. 12 we plot the dispersion 
relation for a 2D model photonic-crystal waveguide and 
compare it with exact numerical results calculated by the 
super-cell plane- wave method [47]. For this case, even 
the simplest tight-binding approximation (i.e., at L = 1) 
gives quite satisfactory results. 

In the tight-binding approximation (L = 1) the dis- 
persion relation can be described by the following simple 
expression 



cos [k(uj)s} — 



D w (uj) ^ lu - u w 
2Vt w (w) uj w A 



where to w is the resonance frequency of the waveguide 
cavities. Furthermore, we have the dimcnsionlcss param- 
eter 



2Vj w (a; w ) _ 

A w — r — ^VlwVvi , 

uj w D' w (lj w ) 



(A7) 



(A6) 



with Vi w = Vi w (uj w ) and v w defined by Eq. (Al), that 
equals half-bandwidth of the waveguide transmission 
band. This band extends from w w (l — A w ) to lu w (1+A w ). 
For our example of photonic crystal waveguide, we find 
A w w 0.052, i.e., its bandwidth is about 10%. 

APPENDIX B: EFFECT OF LONG-RANGE 
INTERACTIONS 

1. Effects of nonlocal dispersion 

As follows from the results of Sec. IIIB above, the 
local-coupling approximation provides us with an ex- 
cellent qualitative analysis of the structure shown in 
Fig. 1(b). However, certain physically important effects 
may be missed in this approximation. A detailed anal- 
ysis of the effects of nonlocal coupling was performed in 
Ref. [22], so that here we may discuss this issue very 
briefly, and may specify it directly to photonic-crystal 
devices. 

In Fig. 5, we provide a comparison of T(uo) calculated 
from Eq. (29) in the local-coupling approximation with 
the exact numerical results for the structure shown in 
Fig. 1(b) for the model photonic crystal described in Ap- 
pendix Al. The results suggest that the local-coupling 
approximation introduces a frequency shift for the band 
edges which agrees well with the corresponding frequency 
shift in the dispersion relation shown in Fig. 12. 

In addition, the resonance frequency is also shifted; it 
is not equal to uj a but is slightly larger. In principle, this 
shift can be produced by two effects: (i) a long-range 
coupling between cavities inside the waveguide and (ii) a 
long-range coupling between the waveguide and the side- 
coupled resonator. First, we explore the former possibil- 
ity. 

Solving Eqs. (24)-(25) for L = 2, we obtain the trans- 
mission and reflection coefficients (l)-(2) with the detun- 
ing parameter 



cr(w) 



Da [ V? w + 3D w V lw V 2w + 3yi w F 2 2 w + 2V 2 3 W cos(3fcs)] - V , a V a , V lvf V2, 
Sm[ J ' Vo, a V a , (VI - VI + D w V 2w ) 



(Bl) 



r 



where all the coefficients are assumed to be frequency- 
dependent analogous to Eqs. (28)- (29). The waveguide 
dispersion k(uS) is now calculated from Eq. (A5) with 
L = 2. 



Fig. 5 shows that the transmission calculated from 
Eq. (Bl) is much closer to the exact numerical results. 
In fact, the nominator of Eq. (Bl) indicates that, indeed, 
the resonance frequency is slightly shifted from the value 
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u) a , and that this shift is proportional to V2w Since V2 W is 
always much smaller than Vi w (see Fig. 11), we can safely 
neglect all the terms proportional to V^w w ith n > 2, and 
obtain the resonance frequency according to 



1 



w 

(V& + 3D w V lw V 2w ) 



(B2) 



Here, the values of all coefhcients are calculated at the 
frequency to a , and v a is dehned by Eq. (Al). 

In addition to the shift of the resonance frequency 
a perfect transmission may occur at the frequencies for 
which the denominator in Eq. (Bl) vanishes: 



2w 



lw 



cos(fcs) ± + cos 2 (fcs) — 2 cos(2fcs) 
1 - 2cos(2/cs) 



(B3) 



However, an analysis reveal that Eq. (B3) has solutions 
only when |1^2w(w)| exceeds |Vi w (w)|/2, a condition that 



appears to be impossible to realize in realistic photonic 
crystals. 



2. Effects of nonlocal coupling 

Another possible reason for a shift of the resonance 
frequency is a nonlocal coupling between the waveguide 
cavities and the side-coupled resonator a. Here, we dis- 
cuss this effect in the framework of the tight-binding ap- 
proximation for the waveguide dispersion (i.e., L = 1) 
to distinguish it from the other type of nonlocal effects 
discussed in the previous sub-section. We assume that 
Vj. a (uj) — V a .j (uj) = for all j > 2, and take into account 
that, for the symmetric structure shown in Fig. 1(b), 
the coupling coefficients are: = Vi. a (uj) and 

V a .-\(uo) = V a .i{uo). Then, we obtain a solution of 
Eqs. (24)-(25) in the form of Eqs. (l)-(2) with the de- 
tuning parameter defined as 



<t(oj) = 2 sin[fcs] 



A»Viw + Vo, a V a ,i + V a ,gVi, a + 2V a ,iVi, a cos(fca) 
[V at o + 2V atl cas(k8)][Vo, a + 2V 1 , a cos(ka)] 

I 



(B4) 



Here, all coefficients are assumed to be frequency- 
dependent, similar to Eqs. (28)-(29). Furthermore, 
the waveguide dispersion k(ui) is calculated again from 
Eq. (A6). 

Eq. (B4) suggests that in this case the resonance 
frequency becomes slightly shifted from the value u> a , 
and this shift is proportional to the values of Vi jCt and 
V a .i, which for our example (see Fig. 5) are equal to 
Va tl = -0.0026 and V lia = -0.0022. Assuming that 
these coupling coefficients are always much smaller than 
V . a and V afi (cf. V afi = 0.096 and V 0<a = 0.082), we 
obtain for the resonance frequency 



L0 o 



, Vo a V a ! + VuoVia 

1 v Vc 

Vlw 



(B5) 



Here, the coefficients are calculated at the resonance fre- 
quency cu a , and v a is defined by Eq. (Al). For the exam- 
ple shown in Fig. 5, this frequency shift is much smaller 
than that described by Eq. (B2) because in this case the 
values of V a< ! and Vi )Q are 3.3 times smaller than the 
value of V2w ■ 



Due to this long-range coupling, there appears a pos- 
sibility of perfect light transmission, as discussed in 
Rcf. [22], but only in the case when |V^ i: l(u;)| exceeds 
|K»,o(w)|/2 or |Vi ia (a;)| exceeds |Vb, a (w)|/2. Again, such 
a scenario appears to be impossible to realize in realistic 
photonic crystal structures. 
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